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ABSTRACT 

We present results of an extensive morphological, spectroscopic, and photometric study 
of the galaxy population of MACS J0025.4— 1225 (z = 0.586), a major cluster merger 
with clear segregation of dark and luminous matter, to examine the impact of mergers 
on galaxy evolution. Based on 436 galaxy spectra obtained with Keck DEIMOS, we 
identified 212 cluster members within 4 Mpc of the cluster centre, and classified them 
using three spectroscopic types; we find 111 absorption-line, 90 emission-line (includ- 
ing 23 e(a) and 11 e(b)), and 6 E+A galaxies. The fraction of absorption(emission)-line 
galaxies is a monotonically increasing(decreasing) function of both projected galaxy 
density and radial distance to the cluster center. More importantly, the 6 observed 
E+A cluster members are all located between the dark-matter peaks of the cluster 
and within ~ 0.3Mpc radius of the X-ray flux peak, unlike the E+A galaxies in other 
intermediate-redshift clusters which are usually found to avoid the core region. In ad- 
dition, we use Hubble Space Telescope imaging to classify cluster members according 
to morphological type. We find the global fraction of spiral and lenticular galaxies 
in MACS J0025 to be among the highest observed to date in clusters at z> 0.5. The 
observed E+A galaxies are found to be of lenticular type with Sersic indices of ~ 2, 
boosting the local fraction of SO to 70 per cent between the dark-matter peaks. Comb- 
ing the results of our analysis of the spatial distribution, morphology, and spectroscopic 
features of the galaxy population, we propose that the starburst phase of these E+A 
galaxies was both initiated and terminated during the first core-passage about 0.5- 
lGyr ago, and that their morphology has already been transformed into SO due to 
ram pressure and/or tidal forces near the cluster core. By contrast, ongoing starbursts 
are observed predominantly in infalling galaxies, and thus appears to be unrelated to 
the cluster merger. 

Key words: Galaxies: evolution, Galaxies:clusters:individual:MACS J0025.4— 1225, 
Galaxies:elliptical and lenticular, Galaxies:starburst, methods:observational 



1 INTRODUCTION 

It is well known that the galaxy population of clusters is 
domi nated by pass i vely evolving , red, early-type galaxies 
(e.g. iDressleil Il980l ; iBower et all 1 1993 ; Ivan Dokkum et al.1 
Il998h . unlike the population of galaxies in less dense 
environments (e.g. iGomez et all 120031 ; [l cwis et all 120021; 
ICooper et al.1 120071 ). Many cluster-specific environmental 
processes, such as ram -pressure strip ping (|Gunn fc Gottl 
Il972l ; iBekki et al.ll2002T ). s trangulation (lLarson et al.lll980l ). 



19961). or tidal dis- 



galaxy-g alaxy harassm e nt dMoore et al.1. 
ruption (|MerrittJ Il983l ; iBvrd fc Valtonenl Il990l ) have been 
proposed to explain the effects of cluster assembly on the 
galaxy population. However, the picture of galaxies being 



accreted in a static cluster environment is overly simplis- 
tic. Many optical and X-ray studies (e.g. Ijones fc Formar] 
1 19991 ; |Pe Propris et al.ll2004) reveal that a significant frac- 
tion, if not all, of the clusters even at recent epochs are 
still growing through mergers of galaxy groups and sub- 
clusters. To get a correct picture of the effects of environ- 
ment on galaxy evolution, the dynamic processes govern- 
ing the formation of these large-scale structures must be 
taken into account. For instance, the predominance of pas- 
sively evolving galaxies in clusters was suggested to have 
its roots in the preprocessing of galaxies in the group en- 
vironment (e.g. |Zabludoff fc M ulchacv 1998; lKodama et al.l 
l200ll ; iTran et al.ll2009l ; IWilman et al.H2009h . thus not being 
directly related to the much denser environment in clusters. 
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Although environment has thus been identified to play 
a critical role in the evolution of galaxies, the effect of events 
that cause the most dramatic change in environment, clus- 
ter m ergers, is still poorly understood (|Girardi fc Bivianol 
2002). For instance, it is still controversial whether star- 
formation activity would be enhanced during c luster merg- 
ers, as favored by numerical simulations (e.g. Bek ki et al.l 
120021 ; [Kronberger et al.l l2008), Obser vationally, this picture 
is s upported by the classical studies of lCaldwell et all(ll993l ) 
and ICaldwell fc Rose! (|l997T ), which find evidence of recent 
star b ursts in merging clusters. More recently, lOwen et al.l 
(2005) suggested that the excess of radio-loud galaxies 
around the cluster Abell 2125 may b e relat ed to the on- 
going cluster merger, and lFerrari et all (|2005T l report on ac- 
tively star-forming galaxies in the region between the two 
sub-clusters of the merging clust er Abell 3921 (z — 0.094). 
In addition, iHwang fc Led (|2009h compare the two merging 
clusters Abell 168 and Abell 1750, and attribute the differ- 
ent amount of galaxy activity, including star formation and 
nuclear activity, in the region between the subclusters to 
the differences in the merger history of the two clusters. On 
the other hand, the ages of numerous post-starburst galaxies 
found in Abell 3921 are too old for th e starbursts to hav e 
been triggered by the ongoing merger (|Ferrari et al.|[2005h . 
Also, a 24^m obse rva tion of C1J0152 7-135 7 (z ~ 0.83, 
lEbeling et aHl2000l ') bv lMarcillac et alj |2007l ) suggests that 
most of the mid-infrared luminous galaxies are infalling and 
thus not related to the merging process. 

The cause of the conflicting observational evidence may 
lie in the complexity of the dynamics of cluster mergers. 
For instance, the specific phase in which a merger is ob- 
served, the viewing angle, and the relative masses of the 
subclusters can all be expected to have a direct effect on the 
triggering of star formation and the resulting spat i al an d 
spectral observational data (see e.g. IHwang fc Led [2009). 
Finding merging systems whose geometry and merger his- 
tory are immediately evident is thus of paramount im- 
portance. To date, only three such systems have been 
identified through a clear segregation between dark mat- 
ter and in tra-cluster gas: the Bul let Cluster (1E0657-56, 
z = 0.296; IClowe et alJliooi l2006t l. Abell 520 (z = 0.201; 
iMahdavi et al.ll2007l ), a nd MACS J0025.4- 1225 (z = 0.586, 
hereafter MACSJ0025; IBradac et all \200i ). Among these 
three systems, MACSJ0025 stands out as the best target 
for a study of the effects of a major cluster merger because 
the masses of the two sub-clusters are similaiQ. 

In this paper, we will report spectroscopic and morpho- 
logical results for galaxies around MAC SJ0025, a cluster dis- 
cover ed in the MAssive Cluster Survey (jEbeling et al.ll200ll . 
120071 ). and, at L x ,chandra = 8.8 ± 0.2 x 10 44 e rgs s" 1 , one of 
the m ost X-ray luminous clusters at z> 0.5 (jEbeling et al.l 
120071 '). Our paper is structured as follows. [J2] summa- 
rizes the dynamical structure of MACSJ0025 as studied by 
IBradac et af] (|2008l V £2 gives an overview of the observation 
and data reduction procedures. ^describes the analysis of 
photometric, spectroscopic and morphological data. In 



1 Note, however, that the supersonic head-on collision of two sys- 
tems of very different mass in the Bullet Cluster enabled a detailed 
study of the effects of ram pressure fr om supersonic gas during a 
cluster merger (see IChung et al1l2009l 'l . 



we present our results. Finally, we discuss a simple merger 
model in ij6] and provide a brief summary in <J7] Through- 
out this paper, we adopt the concordance ACDM cosmology 
with ho — 0.7, Qa = 0.7, Q m = 0.3. Magnitudes are quoted 
in the AB system. 



2 SUMMARY OF DYNAMICS OF MACSJ0025 

MACSJ0025 is the second cluster that shows a significant 
segregation between the hot gas traced by the X-ray surface 
brightness, and the mass traced by the gr avitational shear 
|Kartaltepe et a l. 2008; Bradac et al. 2008) in a simple two- 
body merger configuration. We here briefly summarize th e 
dynamics of the system as presented bv lBradac" et al.l(|200S| y 

Like many other merging clusters, MACSJ0025 features 
a galaxy distribution in redshift space as well as in projec- 
tion on the sky which, taken alone, does not show obvi- 
ous evidence of substructure. Its galaxy redshift distribu- 
tion is consistent with a Gaussian centered at z c ; = 0.5857 
with a c i = 835^59 km s _1 with only mild evidence of a 
bim odal distributiorjfl . In addit ion, the Dressler-Shectman 
test l|Dressler fc Shectmanlll98^ ) detects only marginal clus- 
ter substructure. The redshift difference between the two 
sub-clusters is small (Az = 0.0005 ± 0.0004), which implies 
that the two subclusters are colliding almost in the plane of 
the sky (within 5°). 

The X-ray temperature of the intracluster medium 
(ICM) within a circular region of radius 660 kpc centered on 
the X-ray peak is Q.2&t° Ti keV; the metallicity of the ICM is 
Z = 0.37 ± 0.10 solar. The corresponding X-ray mass of the 
gas is 0.55 ± 0.06 x 10 14 Mq. The virial radius for a cluster 
with the same gas temperature is then ~ 1.3 Mpc, calculate d 
following the prescription provided bv lArnaud et alj {2002). 
Finally, the radius at which a Milky-way-like galaxy would 
lose all its gas due t o ram-pressure stripping is ~ 0.7 Mpc 
(see lMa et all (|2008f ) for details). 

The gravitational lensing analysis of IBradac et al.l 

(2008) shows the two peaks of the mass distribution be- 
ing offset from the X-ray peak by 0.82' for the southeastern 
peak and 0.50' for the northwestern peak. Since the X-ray 
centre lies close to the projected line connecting the gravita- 
tional lensing peaks, the separation of the two lensing peaks 
is thus 1.32' (0.52 Mpc). The masses estimated from the 
lensing model are 2.5±\-° 7 x 10 14 M and 2.6t°A x M 
for the southeastern and northwestern peak, respectively. 



3 OBSERVATION AND DATA REDUCTION 
3.1 Imaging data 

3.1.1 Groundbased observations 

Our object catalogue is based on images in five optical 
bands (B, V, R c , I e , and z') obtained with the Suprime- 
Cam (Mivazak i et al.ll2002h wide-field imager on the Sub- 
aru 8m telescope on Mauna Kea. We also include imag- 
ing data for a near-ultraviolet band (u*) obtained with the 
MegaPrime camera, and two near-Infrared bands (J and 

2 In i|3.2.2l we discuss updates to the spectroscopic sample; the 
results are similar. 
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Figure X. Spatial distribution of galaxies observed spectroscopically: small dots mark the location of all sources with measured spec- 
trum, large dots m ark confirmed cluster members. The overlaid contours show the projected galaxy density as computed by asmooth ( 
lEbeling et alj|2006l . see i|4,2l for details). The red boundary indicates the outline of all DEIMOS masks and defines the study region for 
this paper. The blue dash square indicate the HST ACS image mentioned in £|3.1.2| and defin es the area for the morphology analysis in 
£14.41 The two plus signs show the location of gravitational lensing peaks jBradac et alj|2008l) , and the cross sign shows the location of 
the X-ray flux peak. 



K s ) obtained with the WIRCAM l|Puget et af]|2004 ) instru- 
ment, both on the CFHT 3.6m telescop e on Mau n a Kea . 
We use the sam e pipe line described in iMa et all (2008); 
iKartaltepe et all (120081 ) to reduce the optical and u*-band 
data. Near-infrared images are obtained by stacking t he pre- 
processed data provide d by CFHT using software |Bertinl 
l2006l ; iBertin et al.ll2002r i developed by E. Bertirfl The con- 
struction of a photo metric obj e ct c atalogue then follows the 
same procedure as IMa et all (|2008l ) and IKartaltepe et all 
(2008). The photometry of the Suprime-Cam imaging data 
is calibrated using hundreds of stars with magnitudes rang- 
ing from 16 mag to 19 mag in a nearby SDSS fielcQ. The cal- 
ibration of our MegaCam data is supplied b y CFHT's Elixir 
Finally, the 2MASS catalogue (|Skrutskie et. all 
is used to calibrate the photometry of the stacked 
WIRCAM images. 



3.1.2 HST imaging 

The central ~ 1.3 x 1.3Mpc 2 region of the cluster (Fig. [1} is 
covered by HST ACS images in the F555W filter (4140s) and 
in the F814W filter (4200s), obtained as part of the Cycle- 
13 proposal GO-10703 (PI Ebeling). The data were reduced 
using the HAGGLeS pipeline (Marshall, P. et al. 2009 in 
preparation), which employs MultiDrizzr^H (Koekemoer et 
al. 2002) for the distortion correction, cosmic ray rejection, 
and stacking. Following manual masking of satellite trails, 
cosmic ray clusters, and scattered light in the flat-fielded 
frames, we iteratively refine relative shifts and rotations by 
cross-correlating the positions of compact sources. The indi- 
vidual frames in each filter were then drizzled together onto 
a common 0'.'03 pixel grid using a square kernel with pixfrac 
0.8, applying updated bad pixel masks and optimal weights 
(Schrabback et al. 2009). 



3 http : //astromatic . iap. fr /software 

4 The SDSS photometry is transformed into the Johnson- 
Cousin system using the equation of Lupton (2005) provided on 
http : //www. sdss . org/dr5/algorithms/sdssUBVRITransf orm.html . 

5 http:/ /www. cfht.hawaii.edu/Instruments/Elixir/ 

6 This publication makes use of data products from the Two Mi- 
cron All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis Cen- 
ter/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 



3.2 Spectroscopy 

The spectroscopic data set used here is compiled primarily 
from observations of nine multi-object spectroscopy (MOS) 
masks with the DEIMOS spectrograph on the Keck-II tele- 
scope. The area covered by these masks is marked in Fig. [1] 
We also use 16 additional spectra obtained with the Low 
Resolution Image Spectrometer (LRIS) on the Keck-I tele- 
scope, and the Gemini Multi-Object Spectrograph (GMOS) 
on the Gemini telescope, of objects which meet our selection 
criteria ( i]3.2.1|l and have not been observed with DEIMOS. 



MultiDrizzle version 2.7.0 
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Following the same strategy as iMa et all l|200Sh . 
DEIMOS was used with the 600ZD grating, GC455 fil- 
ter, and centre wavelength at 670(M as a compromise be- 
tween the requirements of wavelength coverage from the 
[Oil] AA3727 line to the Rp line at a redshift of about 0.59, 
and our desire for at least moderately high spectral reso- 
lution. We integrated for 3 x 1800s for each mask; no flux 
calibration was performed. The data were reduced using the 
standard DEIMOS p ipelinjf] developed by the DEEP2 team 
jCooper et al l 120071) ■ Redshifts were determined and veri- 
fied manually using at least two prominent spectral features, 
such as (in absorption) Calcium H and K, H<s , or the G band, 
and (in emission) [Oil] AA3727, Hp, or [OIII] AA4959, 5007. 

3.2.1 Object Selection and Completeness 

In order to maximize the efficiency of cluster member de- 
tections while minimizing potential selection bias, we chose 
our spect r oscop ic targets following the approach taken by 
IMa et all |2008l ). Since, at the beginning of this project, 
imaging data in too few passbands were available to de- 
termine credible photometric redshifts, targets for spectro- 
scopic follow-up observations were selected using the color 
magnitude diagram for the V and R c filters (Fig. [5}. 
To maximize completeness, we limited our survey to rela- 
tively bright objects (itih c ^ 23.0, which is, equivalently, 
Mhc ^M*+2.7 at z=0.59) while adopting a generous color 
cut around the cluster red sequence, extended toward the 
blue end of the distribution. We used our catalogue of pho- 
tometric redshifts to compute the fraction of missed 
blue cluster members(2.3%). We note that incompleteness of 
our spectroscopic survey is severely color-dependent which, 
as we will show in H5.3I could potentially cause the fraction 
of starburst galaxies, but not the fraction of E+A galaxies, 
to be underestimated. 

We define the completeness of our spectroscopic survey 
as the ratio of the number of sources observed to the num- 
ber of sources in the photometric catalogue that meet the 
selection criteria (i.e. that fall within the color band marked 
in Fig.[5]and exceed the magnitude limit), including spectra 
from which we failed to measure a credible redshift. In ad- 
dition, we define the efficiency of our survey as the fraction 
of the number of high-quality spectra from which both the 
redshift and the equivalent width of spectral features can be 
measured accurately (again within the color limits marked 
in Fig. [2}. The completeness and efficiency as a function 
of galaxy magnitude are shown in the top panel of Fig. [3] 
Inside the entire study region, the completeness is roughly 
constant at about 75% at R c < 21.5, with the curve only 
dropping below 40% at the very last bin before the mag- 
nitude limit at R c = 23.0. Also shown in the top panel of 
Fig. [3] is the completeness within a radius of 1 Mpc from 
the centre of the cluster: the completeness is above 80% at 
R c < 22. E0, and drops to 40% at the last bin. Because of 
the significantly reduced completeness at the faint end of 



8 The analysis pipeline used to reduce the DEIMOS data was 
developed at UC Berkeley with support from NSF grant AST- 
0071048. 

9 The completeness may drop because of small number statistics 
at the brightest magnitudes. 



our survey, we adopted a more restrictive magnitude limit 
(R c < 22.5, shown as the dashed lines in Fig. [3]) for the sta- 
tistical analysis of the properties of different spectral types 
of galaxies f i)5.2.2|l to prevent a potential bias caused by 
the luminosity dependence of the spectral type distribution. 
For all other analyses, however, we used the original magni- 
tude limit, since the statistics in the faintest magnitude bin 
(28 cluster members out of 60 galaxies with spectroscopic 
data) are sufficient for our purposes. The differences between 
the completeness curves shown in Fig. [3] (top) demonstrate 
that the completeness of our survey is not spatially uniform; 
the relative spatial distribution of galaxies of different types 
should, however, be unaffected. 

The bottom panel in Fig.[3]shows our efficiency in deter- 
mining the galaxy redshift as well as a secure spectroscopic 
classification. The criteria used to classify galaxies accord- 
ing to spectral type (defined in the following two sections) 
are much more demanding than those required for a simple 
redshift measurement; still, the efficiency for spectroscopic 
classification remains about 0.9 at all magnitudes except for 
the final bin, reflecting the good signal-to-noise ratio of most 
of the spectra. About 10% of our high-quality spectra yield 
spectroscopic redshifts but do not allow a spectral type to be 
assigned; this is usually caused by the key spectral features 
falling on the chip gap or coinciding with telluric lines. At 
the faintest magnitudes probed by this work, the spectral 
classification begins to be limited by the uncertainty of the 
equivalent-width measurements, although we are still able 
to classify about 70% of the respective subsample. 

3.2.2 Spectroscopic catalogue and galaxy redshift 
distribution 

Much of the spectroscopic galaxy sam ple compiled for this 
work was used before in the study of IBradac et alj (|200&t ) 
to perform a basic dynamical analysis. Since then, we have 
doubled the number of galaxies with spectroscopic data; the 
resulting final spectroscopic catalogue comprises 436 galax- 
ies, of which 208 are cluster members (Table[T}. Fig.[T]shows 
the distribution of the spectroscopically observed galax- 
ies within our study region which is highlighted by the 
outline of the DEIMOS masks. The full redshift distribu- 
tion, shown i n Fig. [4] is consiste nt with the one presented 
in Fig. 1 of IBradac et al.1 (|2008l ) and is well described a 
single Gaussian. In addition, the systemic cluster redshift 
(z c ; = 0.5857) and the velocity dispersion (a c i = 843 km 
s _1 ) of the 116 cluster members within the virial radius 
are al so co nsistent with the valu es reported in IBradac et al.l 
(|2008h and lEbeling et all (|2007h . We therefore assign clus- 
ter membership to all galaxies with 0.574 < z < 0.597, the 
redshift range given by z — z c i ± 3<r. 



4 ANALYSIS TECHNIQUE 
4.1 Photometric redshift 

To derive photometric redshifts for the galaxies with 
m Rc <24.2 (~ M*+4.0 at z=0.59V we use the adap- 
tive SED-fitting code Le Phar^\ (|Arnouts et al.1 Il999l ; 

10 http : //www. cf ht .hawaii . edu/~ arnouts/LEPHARE/cf ht_lephare/lephare .html 
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Figure 2. Color magnitude diagram for galaxies inside our study region as defined in Fig[T] Blue filled circles mark objects with observed 
spectra; red filled circles denote spectroscopically confirmed cluster members (cf. Fig. U). Black dots are all remaining objects in our 
SExtractor galaxy catalogue. The dashed lines indicate the selection criteria used for our spectroscopic survey. The red sequence marked 
by the solid line is obtained from a linear fit to the data points redder than V-R c ~ 0.9 and brighter than R c ~ 23.0. 



Table 1. Number statistics of spectral types 



Type 




All 


cluster members 




R c < 23.0 


R c < 22.5 


Rc < 23.0 


Rc < 22.5 


redshift mea- 


436 


383 


212 


182 


surements 










spectral type 


421 


375 


207 


180 


classification 










Emission-line 


236 


204 


90 


75 


galaxies 










Absorption- 


171 


160 


111 


102 


line galaxies 










E+A galaxies 


14 


11 


6 


3 


e(a) 


57 


53 


23 


22 


e(b) 


38 


26 


11 


6 


e(c) 


141 


125 


56 


47 



lllbert et aUbood . |2009). In addition to its \ 2 optimization 
during SED fitting, Le Phare adaptively adjusts the photo- 
metric zero points by using the sample with spectroscopic 
redshifts as a training sell 11 !, which effectively reduces the 
fraction of catastrophic errors and also removed some sys- 




Figure 3. Completeness (top) and efficiency (bottom) of our 
spectroscopic survey as a function of galaxy magnitude. The color 
criterion of Fig. [2]has been applied. Top: survey completeness for 
the entire study region (filled symbols) and for the cluster cen- 
tre region with radius 1 Mpc (open symbols). Bottom: survey 
efficiency, i.e. success rate for redshift determination (open sym- 
bols) and spectral classification (filled symbols). In both panels, 
the dashed line indicates the magnitude limit of the spectroscopic 
sample. 



The resulting adjustments, which reflect the systematic un- 
certainty of the absolute photometric calibration as well as the 
complexity of matching SED templates to the observed spec- 
tra, for all bands are u*:0.081, B:0.066, V:0.090, R c :-0.016, 



I c :0.025, z':-0.0 91, J:-0.181 and K s :0.102, comparable to the 
values quoted by lllbert et al.l (I2009T) . 
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Table 2. Definition of spectral types 



Type 


Criteria 


Emission-line 


detection of [Oil] or Hp ern 


Absorption-line 


no detection of [Oil] and Hp tem , and 




EW(Hi)<4A 


E+A 


no detection of [Oil] and Hp em , and 




EW(H 4 )>4A 


e(a) 


Emission-line galaxies with EW(l/j)> 4A 


e(b) 


Emission-line galaxies with EW(ffj)< 4A 




and EW([OII])< -25A 


e(c) 


Emission-line galaxies with EW(Hg)< 4A 




and EW([OII])> -25A 



Figure 4. Rcdshift distribution for all spectroscopically observed 
galaxies. The systemic redshift of MACSJ0025 is found to be 2 c j = 
0.5857. The dashed lines at z = 0.574 and z = 0.597 (z = z i±3ct) 
mark the spectroscopic extent of the cluster. A zoomed-out view 
of the full redshift range of the galaxies observed by us is shown 
in the sub-panel. 




Figure 5. Photometric versus spectroscopic redshifts for galaxies 
in the field of MACSJ0025. The lower panel shows the residuals 
dz=z p hot _ z a pec as a function of the spectroscopic redshift, while 
the insert in the top left corner shows the histogram of the resid- 
uals as well as a Gaussian model from which we derive Az=0.019 
as an estimate of the uncertainty of the photometric redshifts. 



tained from our database of spectroscopic redshifts by select- 
ing all galaxies with an R c -band magnitude between 20 and 
22, and a minimal separation from the closest SExtractor- 
detected neighbour of 2" to avoid photometric errors from 
close pairs or blended sources. The resulting training set 
comprises 290 galaxies at redshifts ranging from 0.2 to 0.9. 
The comparison of spectroscopic and photometric redshifts 
is shown in Fig. [5] The statistical error of the photomet- 
ric redshifts, as derived from the Gaussian distribution of 
the residuals shown in the insert of Fig. [S] is found to be 
Az = 0.019. 

Based on the photometric redshifts thus computed, 
we define cluster members to be all galaxies with 
\z ph - z c i\ < 2a ph -z, where a ph - z = (1 + z c i)Az. Note that 
this redshift cut (±2a p h-z = ±0.061) is much more gener- 
ous than the one used before for the spectroscopic redshifts 
(±3<t = ±0.012); the chosen cutoff value represents a good 
compromise between completeness and contamination. The 
resulting catalogue of "photometric" cluster members in- 
cludes 1602 galaxies within the entire region imaged with 
SuprimeCam and allows us to estimate the local galaxy den- 
sity within and outside our spectroscopic study area. 



4.2 Projected galaxy density 

Using the photometrically selected set of cluster members, 
we calculate the projected galaxy density in two ways. For 
the convenience of visual demonstration, a smoothed ver- 
sion of the spatial distribution of cluster members (shown 
in Fig.[T]) is obtained with the adaptive-smoothing algorithm 
asmooth IjEbeling et al.l l2006) which preserves the signal-to- 
noise ratio of all significant features on all scales. In ad- 
dition, in order to allow a direct comparison with the re- 
sults from earlier w ork on the impact of cluster environment 
on galaxy activity dKodama et al.ll200ll ; iBalogh et al. 1 12004 
Sato fc Martinll2006l ; iMa et al.ll200cf ), we also calculate, for 
all "photometric" cluster members, the projected local den- 
sity Eio , denned as the galaxy number density within a circle 
whose radius is given by the separation from the 10 th closest 
neighbour. The densities obtained from the two approaches 
are consistent with each other. 



tematic trends in the diffe rences between th e spectroscopic 
and photometric redshifts (|llbert e t al. 2006). We us e the li- 
brary of empirical templates from Ilbert et alj (|2009T ). which 
is included in the Le Phare package. The training set is ob- 



4.3 Spectroscopic analysis 

We perform equivalent- width (EW) measurements and spec- 
tral type classification following the prescription given 
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Figure 6. ACS images (F505W and F814W) of selected cluster galaxies with stamp size 4.0". The six E+A cluster members are shown 
in the top two rows, and randomly selected emission-line and absorption-line galaxies in the cluster are shown in the middle and bottom 
two rows respectively. Note that the first E+A galaxies on the top row may be contaminated by a blue lensing arc. 
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in iMa et all f2008j3- I n short, we classify galaxies into 
three major spectral types (absorption-line, emission-line, 
and E+A galax ies), and three sub- types (e(a), e(b), and 
e(c)) (see also jPressler et alj [l999l: IPoggianti et all 1 19991 ; 



IPoggianti fc Wvj|2000l ; iTran et alj|2003l ; lGotoll2004 for the 
properties of different spectral types). Specifically, E+A, 
also called "post-starburst" , galaxies are identified by the 
strong Balmer absorption lines indicating significant amount 
of young stars, and the absence of emission-lines indicat- 
ing current star-formation activity has been terminated. On 
the other hand, e(a), also called "dusty starburst", galax- 
ies are mostly useful to trace the starburst galaxies. Their 
strong Hs absorption lines indicate the recent burst of many 
young stellar population, and, in contrast to the E+A galax- 
ies, the moderate [Oil] emission exhibits the starburst has 
not been terminated. Furthermore, the infrared study of 
IPoggianti fc Wul (|2000l ) finds about half of the very lumi- 
nous infrared galaxies [log(L/_R/LQ)>11.5] in their sample 
to exhibit e(a) spectra, which suggests that their moderate 
[Oil] emission is heavily extinguished by dust. Also, the lack 
of e(b) spectra in the same sample indicates that e(b) char- 
acteristics are less well suited to identify starburst galaxies, 
although they exhibit stronger [Oil] emission. For clarity, 
the definition of the spectral types is summarized in Tabled 
The number of galaxies of each type is given in Table [1] 
Note that the fraction of E+A galaxies in this cluster is 
~ 2.9 ± 1.2% for m Rc <23.0, and ~ 2.2 ± 1.1% for the more 
stringent magnitude limit of mR c <22.5. The fraction is rel- 
ativ ely low compared to other intermediate-redshift cluster s 
sec IMa et al. 2008; IPoggianti et al. 2009]; ITran et al.ll20ol ) . 



4.4 Morphological analysis 

Taking advantage of the availability of high-resolution imag- 
ing data obtained with HST/ACS for the central 1 Mpc re- 
gion of MACSJ0025 (the sub-panel in Fig.[T|, we attempt to 
examine the morphological transformation of galaxies dur- 
ing the cluster merger. Within this region, morphological 
types are assigned to all cluster members with mR c < 23.0, 
including 80 spectroscopically identified galaxies, as well as 
28 galaxies classified as cluster members based on their pho- 
tomotric redshifts, using three broad categories, namely el- 
liptical (E), SO, and late-type (S). The classification is per- 
formed visually by examining galaxy images in the F814W 
band (in random orde r) using templa t es fro m similar clas- 
sification efforts, e.g., IPostman et alj (|2005l ). Fig. [6] shows 
postage stamps for a selected subset of galaxies with the 
visual morphological classification noted. To test the ro- 
bustness of the visual classification and to obtain a quan- 
titative measure of galaxy morphology, we also study the 
bumpiness-Sersic n distr ibution for our galaxie s, following 
the technique outlined in Blakcslc e et alj (|2006l ). 



12 The only difference is that we use the Hj abso rption line, 
rather than both Hg and H 7 as in IMa et alj feOOSt) . to classify 
E+A galaxies, since, for galaxies at the redshift of MACSJ0025, 
H 7 falls onto the strong O2 absorption feature of the atmosphere. 



4-4-1 Surface-brightness fitting 

We fit the two-dimensional surface-brightness distribution, 
as obtained from ACS F814W images of our galaxies, with 
a Sersic model and a n additional sky- background compo- 
nent using GALFnPI (|Peng et alj|2002l) . All parameters of 
the Sersic model, including the coordinates of the galaxy 
centre, normalization, Sersic index, effective radius, ellip- 
ticity, and orientation are fitted simultaneously. Following 
IBlakeslee et~ai1 (|2006h . we constrain the Sersic index, n, to 
0.2 < n < 4.0 since larger values do not substantially im- 
prove the fit but sometimes cause the best-fit value of the 
effective radius to become unreasonable large. For the initial 
values of the fit parameters, critical for GALFIT, we adopt 
numbers derived using SExtractor, except for the Sersic in- 
dex, which cannot be calculated directly using SExtractor 
and is set to 2.5 for all galaxies. We compare and examine 
all fit results, including the residuals, to ensure that the al- 
gorithm did not get trapped in a local minimum. If the fit 
does not converge or produces implausible results, the initial 
parameters are adjusted manually. 

GALFIT is run on images whose size is determined by 
8 FLUX_RADIUS estimated using SExtractor, but requiring 
a minimal size of 5". A large image size is advatageous in 
that it provides better estimates of the sky background and 
the galaxy outskirts. On the other hand, larger images will 
also include more neighbouring objects, especially near the 
dense cluster core. We address this issue by masking out 
any close-by objects which are not connected to the target 
galaxy in the segmentation map generated by SExtractor. 
If a second galaxy is connected to the target galaxy in the 
section map, it is fitted with a separate Sersic model. 

Galaxies visually classified as irregular are excluded 
from Sersic model fits. 



4-4-2 Bumpiness 

The Sersic index is a well-defined parameter quantifying the 
concentration of the bulge; it correlates with morphology 
such t hat early-type galaxies usually have larger Sersic in- 
dices l|Andredakis et ai1 IT995). However, the range of Sersic 
indices observed for a given morphological type is too large 
for it to be reliable morphological indicator in its own right. 
One of the reason for which Sersic indices cannot be mapped 
accurately to morphologial types is that, in reality, galaxies 
are much more complex than the smooth Sersic model: more 
than one Sersic model may be needed to describe the bulge 
and the disk, and non-radial features, such as spiral arms or 
tidal tails (crucial for visual classifications) cannot be mod- 
elled at all. To overcome t his limitation, various alternative 

T99l ; iTakamival 



and additional techniques llAbraham et al 
Il999l; IConselicd 120031 ; IBlakeslee et alj |2006) have been in- 



troduced. From these a pproaches, we have cho sen the B—n 
classification defined in IBlakeslee et~aT1 i|200rj) . because of 
its promising performance in distinguishing SO galaxies from 

elliptical and spiral galaxies. 

Following Blake slee et al] (|2006l ). we define the bumpi- 
ness (B) as the normalized residual of the surface-brightness 
distribution, I, and the fitted Sersic model, S(R e ,n): 



13 http : //users . obs . carnegiescience . edu/peng/work/galf it/galf it .html 
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Figure 7. Distribution of Sersic indices and bumpiness parameters of MACSJ0025 cluster members. The symbols indicate the mor- 
phological classes determined visually: early-type (circles), lenticular (SO, triangles), and late-type galaxies (stars). The symbol size is 
proportional to the total brightness in R c . The two solid lines show the linear rel ations B = 0.065(re + 0.8) and B = 0.05(n — 1), which 
divide early- and late- type galaxies, and SO and elliptical galaxies, respectively {Blakcslce et al. 2006). In the left panel, the shading 
indicates morphological sub-type (black represents E, very dark grey E/SO, dark grey SO/E, gray SO, and so on). Filled symbols identify 
cluster members with spectroscopic redshifts; open symbols mark galaxies classified as cluster members based on their photometric 
redshift. In the right panel, we show only cluster members that have spectroscopic data; here the color indicates the spectroscopic types: 
blue for emission-line, green for E+A, and red for absorption-line galaxies (see Section 15.1.11 . 



B = 



y/([/-g(fl e ,n)5-(af) 
(S(R e ,n)) 



(1) 



5 RESULTS 



5.1 Morphology 



where R e and n are the effective radius and index of the 
Sersic model, and a is the flux uncertainty of the observed 
surface-brightness image. We compute these averages over 
an annulus with an inner radius of two pixels from the galaxy 
centre and an outer radius of 2R e . The subscript s indicates 
smoothing by a Gaussian function with FWHM= 0.085". 

To compare the "visual" morphology and the B—n clas- 
sification, we plot bumpiness against Sersic index for our 
morphological sa mple (left panel of Fig . [3). We find that the 
division lines of iBlakeslee et all (|2006 ) are also applicable 
to our sample, supporting the notion of Ivan der Wei et al.l 
(|2007h that the B—n technique is robust and does not re- 
quire adjustments for different studies. Several outliers are 
notable in Fig [7] Individual inspection reveals that, for in- 
stance, two galaxies visually classified as early type (red cir- 
cles near the upper division line) are both Brightest Cluster 
Galaxies (BCG s), which cannot b e fitted well by a po wer-law 
profile (see e.g. lLaine et all 120031 : lLauer et alj|2007t > result- 
ing in B—n parameters close to that of late-type galaxy. On 
the other hand, the isolated galaxy close to the lower division 
line (bumpiness B ~ 0.12 and Sersic index n ~ 3.4), visually 
identified as a face-on late-type spiral with faint but obvious 
spiral arms, exhibits a very faint disk, but a very prominent 
bulge, a combination which results in a small B value. 

Indeed. felakeslee et alj l|2006l ) note that the B—n classi- 
fication scheme does not always identify SO galaxies reliably, 
and that the agreement with visually identified SO is about 
60%. This is consistent with our own findings: about 60% of 
SOs (as identified by their location in the B—n diagram) are 
also visual SOs, and about 80% of visual SOa fall between 
the dividing lines of Fig. [7] We will, in the following, use 
both classifications to check the impact of this uncertainty 
on our conclusions. 



5.1.1 Correlation between morphology and spectral type 

We compare the morphology of cluster members of different 
spectral types in the right-hand panel of Fig. [7J Not surpris- 
ingly, almost all of the "visual" late-type galaxies feature 
emission lines whereas almost all of the "visual" SO and 
early-type galaxies show pure absorption spectra. Notable 
ex ceptions are three so-called passive spi r als (first found 
in ICouch et all 1 19981 : iDressler et ail 1 19991 : iPoggianti et all 
I1999T I and, conversely, one elliptical and three lenticular- 
galaxies (again, visually classified) which exhibit emission 
lines in their spectra. 

Of particular interest to us, in the context of this work, 
are the E+A galaxies. We find them to be very similar to 
the population of lenticular galaxies in all respects, i.e., vi- 
sual morphology (see Fig. [BJ, value of the bumpiness pa- 
rameter, and in the relatively narrow range in Sersic index 
values (n ~ 2) that they exhibit. Almost since their di scov- 
ery in clusters (ICouch et al.lll994l : [Dressier et al.lll994T ), the 
morphology of E+A galaxie s has been a topic of particular 
interest. Many studies (see iGotd ( 20051 ) as well as the re- 
cent review in IVergani et al. (201C)) find E+A galaxies to 
be mostly bulge dominated, although evidence of interaction 
can usually be found, especially for post-starburst galaxies 
in the field. The findings presented by us here are consis - 
tent with earlier work (e.g. for Coma. [Poggianti et al.ll2004r ) 
in that we also find E+A galaxies to be bulge-dominated 
systems, although their surface-brightness profiles are sig- 
nificantly flatter than those of other absorption-line galax- 
ies. Visual examination of the post-starburst galaxies in 
MACSJ0025 does, however, not reveal any evidence of re- 
cent interaction, a finding that is supported by the bumpi- 
ness values measured which are much lower than those of 
late-type galaxies. These trends (illustrated in Figs. [7J and 
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Figure 8. The distribution of galaxies of different spectral types in the cluster centre. Symbols are as follows: triangles mark lenticular 
galaxies, stars mark spirals, and circles mark elliptical galaxies. All symbols are also colour-coded, such that blue denotes emission lines 
in the galaxy spectrum, red absorption lines only, green marks E+A galaxies, and white flags galaxies identified as cluster members only 
through their photometric redshifts. The three large circles, finally, delineate the three study regions labelled as discussed in £15.1.21 



[6| add to a picture in which E+A galaxies in clusters exhibit 
different properties, and thus most likely follow different evo- 
lutionary paths, than their counterparts in the field (see also 
iTran et alj|2003 : iGotdliooEl ; IPoggianti et alj|2009h . 



5.1.2 Spatial distribution of morphological types 

We calculate the fraction of morphological types in four dif- 
ferent regions: 1) 1 Mpc (radius) around the X-ray lumi- 
nosity peak, i.e. essentially the entire ACS field of view, 2) 
0.15 Mpc around the X-ray luminosity peak, 3) 0.15 Mpc 
around the peaks of the lensing-derived mass distribution, 
and 4) region (1) again but this time excluding regions (2) 
and (3) (see Fig. |SJ. The results are summarized in Ta- 
ble [3] and Fig. [9] and are consistent for our visual and 
the B — n classification. We note the following trends: a) 
elliptical galaxies dominate the regions around the peaks 
of the mass distribution, as expected from the well known 
morphology-density relation l|Dressler|[i980h : b) the region 



between these peaks, however, features an elevated fraction 
of SOs, 0.67 ± 0.14, inconsistent with the global value for 
region (1) at the 90% confidence level. 

In general, the global fractions of the three morpho- 
logical galaxy types in MACSJ0025 fall within the range 
of value observed in other clusters at similar redshifl (e.g . 
Postman et al.ll2005l:|Desai et al.ll2007l ; iBlakeslee et~aHl2006l : 



Poggianti et al . 2009). To facilitate a comparison with pre- 



vious work we recompute our results for a region defined by 
0.6 R vir = 0.75 Mpc (for MACSJ0025) and a magnitude 
limit (M„ < —2 0.0), thereby adopting the constraints used 
in the study by IPoggianti et all (|2009l i. We find our results 
for MACSJ0025 to, roughly, follow the trends with velocity 
dispersion X-ray luminosity observed for other distant clus- 
ters (see Fig. I10|) . but note that systematic effects are likely 
to add to the already substantial statistical uncertainty: in- 
fall along the line of sight (a well known bias for optically 
selected clusters) will temporarily inflate the velocity disper- 
sion, and merger activity along any axis will boost the X-ray 
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Figure 9. Fraction of early-type and lenticular clusters members 
in different regions of the clusters. The open symbols show results 
based on visual morphological classification; values marked by 
filled symbols were obtained using the B-n classification scheme. 
The labels identify the four study regions following the convention 
established in Fig.|8]and described in the text. The plotted error 
bars represent la confidence. Note the high fraction of elliptical 
galaxies near the lensing-derived mass peaks (region 3), and (less 
statistically significant) the increased fraction of lenticulars near 
the X-ray peak in between (region 2). 
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Figure 10. Reproduced from [Poggianti et al Black cir- 

cles denote nearby clusters (0.04 <z< 0.07), and red circles mark 
distant clusters (0.5 <z< 1.2). The open (visual classification) 
and filled (B-n classification) blue circles show our morpholog- 
ical fractions for MACSJ0025 evaluated within the same radius 
of 0.6 K v ir- 



Figure 11. Fraction of galaxies of different spectral types as a 
function of projected local density. The symbols denote the distri- 
bution of different spectral types: the red-circles for absorption- 
line galaxies, the green-circles for E+A galaxies, blue-squares 
for emission-line galaxies, cyan-squares for e(a), and magenta- 
diamonds for e(b). Error bars assume Poisson statistics. Some of 
the symbols are slightly shifted horizontally for clarity. 
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Figure 12. Distribution of galaxies of different spectral type with 
cluster-centric distance. The symbols the same as in Fig. If II In 
each bin the fraction is calculated by dividing the number of 
member galaxies of a given spectral type by the total number 
of member galaxies observed spectroscopically. The plotted error 
bars assume Poisson statistics. Some symbols are shifted slightly 
by ±0.02 Mpc for clarity. 



luminosity beyond the virial value. With this major caveat 
in mind, we note nonetheless that the fraction of S and SO in 
MACSJ0025 is the highest among the clusters with z> 0.5. 

5.2 Spatial distribution of spectral galaxy types 

5.2.1 Projected galaxy density 

In Fig.[TT]we plot the fraction of galaxies of different spectral 
type as a function of the projected local galaxy density, Eio, 
estimated using all photometrically selected cluster mem- 
bers. The observed monotonic fall and rise of the fractions 
of emission- and absorption-line galaxies (top panel) with 



increasing galaxy density, Eio, is expected from the relation 
between star-formation rate and local density, as well as 
from the morphology-density relation, established by stud- 
ies of o ther clusters at i nterm e diate redshift, e .g. | Dresslerl 
jl980j) ; iPoggianti et all (|l999l ); iKodama etafl (|200ll ). The 
fraction of absorption-line galaxies in the densest regions is 
about 70%, dropping to about 30% at the lowest density 
studied in this paper. Conversely, the fraction of emission- 
line galaxies increases to ~ 70% at the lowest density probed 
here, which is consi stent with the fraction of emission-lin e 
galaxies in the field (|Gomez et al.ll2003l ; lc~ooper et alj|2007t ). 

In the bottom panel of Fig. 1111 we see that all six 
E+A galaxies identified in MACSJ0025 reside in regions of 
very high galaxy density. The distribution of E+A galax- 
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Table 3. Fraction of morphological types 



Region 


# of galaxies 


S 


visual 
E 


SO 


S 


B—n identified 
E 


SO 


1: entire ACS image (IMpc) 


96 


0.36 ±0.04 


0.34 ± 0.04 


0.29 ±0.04 


0.39 ± 0.04 


0.29 ±0.04 


0.32 ±0.04 


2: X-ray centre (0.15Mpc) 


6 





0.33 ± 0.16 


0.67 ±0.14 


0.17 ±0.14 


0.17 ±0.14 


0.67 ±0.19 


3: Lensing peaks (0.15Mpc) 


22 


0.14 ±0.07 


0.55 ± 0.07 


0.32 ±0.08 


0.14 ±0.07 


0.59 ±0.07 


0.27 ±0.08 


4: Reg 1 but excluding Reg 2,3 


66 


0.44 ± 0.05 


0.29 ± 0.05 


0.27 ±0.05 


0.47 ±0.04 


0.21 ±0.04 


0.32 ±0.05 


5: X-ray centre (0.75Mpc)* 


89 


0.35 ±0.04 


0.36 ± 0.04 


0.29 ±0.04 


0.37 ±0.04 


0.30 ±0.04 


0.33 ±0.04 



Following the definition in lPoggianti et al ] l|2009l l: for galaxies within 0.6 R v ir = 0.75 Mpc and magnitude M v < —20.0. 
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Figure 13. Projected spatial distribution of different galaxy types in the MACSJ0025 system. The symbols are the same in Fig. 1111 In 
addition, black circles mark the locations of the cluster members with spectroscopic redshift and small black dots mark the locations of 
all galaxies within our study region that have been observed spectroscopically. Also shown are the galaxy-density contours from Fig. [T] 
The red circles have radii of R rs , as estimated in [[2] and a scale of 1 Mpc is drawn on the bottom- left corner of the figure. 



ies in MACSJ0025 is thus very different from that observed 
in other intermediate-redshift clusters (e.g. iDressler et alJ 
1 19991 ; iTran et alj|2003l ; iMa et al.ll2009l) where E±A galaxies 
are found in a wide range of projected galaxy densities and 
do not favor the densest environment. 

In contrast to the distribution of E±A galaxies, the frac- 
tion of e(a) galaxies is the lowest in the densest regions of 
the cluster, increases to a maximum at Eio ~ 10 Mpc -2 
(which happens to be the density at which the fractions of 
emission- and absorption-line galaxies are the same), but 
does not rise further at lower densities. This result is consis- 
tent with the consensus that intermediate galaxy densities 
provide the most efficient environment for star formation. 



5.2.2 Cluster- centric distance 

Plotting the distribution of cluster galaxies of different spec- 
tral type as a function of cluster-centric distance (Fig. 1121 a 
two-dimensional version is shown in Fig. 1 13 p , we find a trend 
that is very similar to the one shown in Fig. 111! i.e. a smooth 
monotonic decrease/increase of the absorption/emission-line 
galaxies with increasing distance from the cluster centre. 
Note, however, that the fraction of emission-line galaxies is 



significantly suppressed only within r ~ 1 Mpc, which is ap- 
proximately the size of the virial radiuj^l 

Similar to the distribution as a function of galaxy den- 
sity depicted in Fig. 1111 the distribution of e(a) galaxies 
as a function of cluster-centric distance also peaks at the 
scale (r ~ 1 Mpc) at which the fractions of emission- and 
absorption-line galaxies are roughly equal, with no e(a) 
galaxies being found in the very centre of the cluster. The 
trend can be seen more clearly in the two-dimensional dis- 
tribution shown in the right panel of Fig. [13] This result 
implies the presence of an efficient trigger of starburst ac- 
tivity near the virial radius of the cluster. 

Most interestingly, however, all of Figs. [8] 1121 and [13] 
illustrate clearly that E±A galaxies appear to reside pref- 
erentially in the very centre region of the cluster (note: we 
define the cluster centre to be the location of the peak of 
the X-ray emission) , a result that is in obvious conflict with 
the consensus t hat E±A galaxi es in clusters avoid the very 
cluster core fsee lMa et al.l 20081 . and references therein). As 



14 Since the cluster-centric radius is difficult to define in this 
merging system we adopted, for this analysis, the location of the 
peak of the X-ray surface brightness as the cluster centre. The 
exact choice is not critical though at the coarse binning used 
in Fig. 1121 We confirmed that our results do not change when 
the cluster-centric radius is computed using a more complicated 
elliptical model where the two focal points are given by the two 
peaks of the lensing mass distribution. 
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galaxies, — 19 < M_g < —16 (e.g. iBarazza et al ,l2009h. found 
in local galaxy clusters and groups (see Binggeli et al.lll99(il ; 
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Figure 14. Redshift distribution of different spectral types of 
cluster members. Redshift distribution of the emission-line and 
absorption-line cluster members are shown, respectively, with 
blue and red histograms, overlaying on the redshift distribution 
of the entire spectrum sample (grey). The magenta, cyan, green 
points, mark the redshift of the e(b), e(a) and E+A galaxies. 

we shall discuss in Section [6] this apparent conflict can be 
resolved once the geometry and dynamics of this merger are 
taken into account. 



5.2.3 Redshift distribution 

The redshift distribution for different spectral types of 
cluster members is shown in Fig. 1141 It is narrower for 
absorption-line galaxies than for emission-line galaxies, 
which is expected if the infalling population is dominated by 
emission-line galaxies. More importantly, we note that the 
redshift distributions of E+A, e(a) and e(b) galaxies are not 
significantly different from that of emission- or absorption- 
line galaxies, consistent with the hypothesis that, in a dy- 
namical sense, they are drawn from the same parent popu- 
lation. 



5.3 Color— magnitude distribution of different 
spectral types 

We plot the V— I c colour and R c magnitude of all cluster 
members in Fig. 1151 using different symbols to denote spec- 
tral type. Similar to previous results (e.g. iMa et al] l200ct 
IPoggianti et~"aH 120091 ; fTran et al.ll2003l ), the E+A galaxies 
are found to inhabit the "green valley" , with one exception, 
for which the (groundbased) colour measurement is erro- 
neous due to contamination by a blue gravitational arc, as 
revealed in our ACS image (top left image in Fig. [6]). All 
E+A cluster members are fainter than M_r > —20.7, unlike 
the E+A galaxies in other intermediat e redshift cluster s, in- 
cluding MACSJ0717.5+3745 for which lMa et~"aT1 (|2008h find 
M R = [-23.5,-20.0]. The magnitudes of E+A galaxies in 
MACSJ0025 are thus closer to those of the faint populations 
of E+ A galaxies found to d ominate in local clusters, e.g. 
Com a (IPoggianti et alj|2004l ) and Abell 3921 (jFerrari et al.l 
2005). In the B band, the range in absolute magnitude of 
these E+A galaxies is Mb = [—19.3,-20.2], which is just 
above the limit of the typical magnitude of early-type dwarf 



lMatedll99cj|Haines ; et "al]|2007l ; lBarazza et al1l2009l . and ref- 
erences therein). Note, however, that the B-band magni- 
tudes of the E+A galaxies in MACSJ0025 are boosted by 
the presence of a young stellar population, and that their 
brightness will drop by 3-4 magnitudej^l. as they evolve 
passively from z = 0.59. Moreover, the morphology of the 
E+A galaxies in MACSJ0025 is well matched to that of 
early-type dwarf galaxies which are characterized by a flat- 
ter surface-bright ness profile (Sersic index of about 2; (e.g. 
ISmith e~a l. 2009)) than for typical ellipticals. 

Unlike E+A galaxies, e(a) galaxies exhibit a broad 
range of magnitudes within the "blue cloud". Because of 
their strong Ha absorption features and the moderate EW 
of their emission lines, e(a) gala xies are believed to be possi- 
ble progenitors of E+A galaxies dPoggianti et al.lll999l ). They 
are, in general, brighter than E+A galaxies, which might, 
again, be due to excess stellar light during the starburst 
phase. In this scenario, their magnitude would drop to the 
levels typical of E+A galaxies once the starburst is termi- 
nated. 

Finally, we note that e(b) galaxies are observed predom- 
inantly at the very blue end of the explored colour range. 
We therefore suspect that many cluster members of e(b) 
type are missed in our study due to the colour criterion ap- 
plied for the selection of the spectroscopic sample; hence, 
the fraction of e(b) galaxies should be considered as a lower 
limit. 



6 PHYSICAL INTERPRETATION 
6.1 A simple model 

Because of its simple merger geometry (major equal-mass 
merger, in the plane of the sky, after first core passage), 
MACSJ0025 provides us with an excellent opportunities to 
study the effect of cluster mergers on fundamental galaxy 
properties, specifically their morphology and star formation 
activity. 

Our understanding of the merger history and dynamics 
of MACSJ0025 is depicted schematically in Fig. [161 The key 
phases and physical processes at work can be described as 
follows: 

Phase A: pre-contact; no significant galaxy-gas interaction 
yet; relative galaxy velocities too high for galaxy mergers, 
but increasing probability of harassment events 

Phase B: first contact; increased ICM interaction induces 
shocks; onset of segregation between dark and luminous 
matter; ram-pressure stripping and tidal interactions begin 
as galaxies encounter steep gradients in gravitational po- 
tential and ICM density; morphological tranformation and 
triggered star formation 

Phase C: core passage; significant ICM interaction and 
shock heating; non-collisional components (dark matter, 
galaxies) proceed as gas is left behind; galaxies on the lead- 
ing edge of either cluster now fully stripped of their gas 



15 Estimated using the simple stellar synthesis model from 
iBruzual fc Charlotl J2003T) 
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Figure 15. Color-magnitude diagram for galaxies of different 
spectral type. The the symbols are the same as in Fig. 1111 The 
E+A galaxies (green filled circles) are fainter than mrt c > 22.0 
and are located mostly in the "green valley" below the red- 
sequence marked by the absorption-line galaxies (red open cir- 
cles). Compare to the E+A galaxies, the e(a) galaxies (cyan- 
squares) are brighter and exhibit a broad range of color in "blue 
cloud" . 



content; galaxies at the trailing edge experience Phase B 
effects 

Phase D: present stage; dark matter cores fully separated 
but single ICM component concentrated at global cen- 
tre of mass; gas-galaxy interactions (triggered star forma- 
tion, stripping) continue as trailing galaxy populations pass 
through core. 

Qualitatively, this picture is supported by the re- 
sult of numerical simulation of ram-pressure strippi n g (e.g 



Bckki 



Bekki 



Fuiita fc Naeashimal Il999l : iKronberger et all l200Sl: 
2009 1 ) and tidal inte raction (|Bvrd fc Valtonenlll99ol ; _ 
19991 ; lGnedin|[2003bl lah . Although the environments encoun- 
tered in this process by galaxies at the leading and trail- 
ing edges of the individual clusters are roughly symmetrical, 
galaxies on the far side of the merger experience a hostile en- 
vironment for longer since they encounter strong tidal fields 
(passage through the dark matter core) and ram-pressure 
stripping sequentially; in addition, they pass through the 
more massive and denser post-collision gaseous core. 

We note that the velocities of galaxies relative to their 
respective cluster core are not negligible, and hence their 
observed relative locations may differ substantially from the 
ones prior to collision. Exceptions are, however, the E+A 
and the e(a) populations, both of which will remain rea- 
sonably faithful tracers of environment thanks to the rela- 
tively short duration of the star-formation and quenching 
processes. 

In the following, we will examine whether the features 
observed in the galaxy population of MACSJ0025 support 
the picture described above. 



6.2 Distribution of E+A galaxies 



As discussed in H5.2.1I and N5.2.2I we find all but one of the 
E+A galaxies in MACSJ0025 to be located near the peak 
of the X-ray surface brightness, and thus between the dark- 




Figure 16. A schematic illustration of the collision process. Red 
circles mark the centres of the dark-matter distributions of the 
two clusters, while light blue and black circles represent the distri- 
bution of the ICM and the cluster galaxy population, respectively. 
See text for discussion. 



matter cores of the individual clusters. This is exactly what 
would be expected in the scenario outlined in £16.11 and de- 
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picted in Fig. 1161 The merging process in MACSJ0025 thus 
appears to cause dramatic and rapid galaxy evolution (the 
triggering and termination of startbusts) in an environment 
where it is rarely observed, namely the region of highest gas 
density in the very cluster core. 

We should mention in this context that E+A galaxies 
have in fact been detected in cluster cores before. In all cases, 
however, projection effects seem to play a major role. Specif- 
ically, the group of starburst galaxies found in Abell 851 (see 
lOemler et al.l 12009 . and references there in) forms a group 
of high peculiar velocity falling into the cluster along the 
line of sight. By contrast, no evidence of projection effects 
is found for MACSJ0025 (see ^2~3|) . 

6.3 Morphological transformation 

In the merger scenario discussed above, the galaxy popula- 
tions of both clusters should have been exposed to condi- 
tions conductive to trigger star formation and ram-pressure 
stripping for a large fraction of the duration of the cluster 
collision. As a result, one would expect the morphological 
transformation from late-type to SO galaxies to have been 
particularly efficient in this system in recent times. Indeed, 
we measure a high global SO fraction of about 0.3 (Fig. 
and Table [3]), among the highest found to date in any clus- 
ter at z > 0.5, but still wi thin the observed range of 0.05- 
0.33 (|Poggianti et al.ll2009h . Furthermore, the high global S 
fraction of about 0.39 may reflect that cluster merger help 
the mix of late- type galaxies at outskirt. A more dramatic 
increase in the global SO and S fraction would have been 
surprising, given that the morphological mix of galaxies in 
clusters is constantly driven to a global average as SOs con- 
tinue to evolve into ellipticals, and late-type galaxies are be- 
ing accreted from the field. A snapshot of the merger-driven 
transformation process is, however, provided by our view of 
the core of MACSJ0025 where we find all E+A galaxies to 
be lenticulars, leading to a local SO fraction of about 70%. 

Our hypothesis that interactions with the gravitational 
potential and the intracluster medium are primarily respon- 
sible for the accelerated evolution of galaxies in this massive 
merger is also supported by a comparison with the mor- 
phology of E+A galaxies in the field. Signature of tidal 
destruction or galaxy merging c an be removed effic iently 
and quickly by both ram press ure l|Kapferer et al-lfeoOST ) and 
tidal forces fBekki et al.ll200ll) in the cluster core. 

6.4 Temporal evolution 

We can use estimates of the duration of characteristic phases 
of galaxy evolution derived from spectral diagnostics for a 
consistency check of our simple merger model (Fig. I lfcif) . 

Specifically, ste llar synthesis models 

jBruzual fc Chariot! I2003T I predict how the strength of 
the D„(4000) break and the EW of Hf vary along an 
evolutionary track from the beginning of the starburst to 
several Gyrs after. In Fig. 1171 the data for galaxies of the 
various spectral types observed in MACSJ0025 are overlaid 
on the theoretically expected trackj^l. We note that these 
models can only be used to describe the absorption features 

16 The models assume a delta-function starburst in a 5 Gyr-year- 
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Figure 17. The D n (4000) and EW of H<5 of galaxies in the clus- 
ter; cyan for e(a) galaxies, green for E+A galaxies, magenta for 
absorption-line galaxies in the regions farther from the cluster 
cores along merger axis, and red for absorption-line galaxies in re- 
gion 3 of Fig. [8] Th e four curves represent ste llar-population syn- 
thesis models from lBruzual &: C hariot (2003]), assuming a delta- 
function starburst in a 5 Gyr-year-old galaxy with exponentially 
decaying SFR (r = XGyr). From top to bottom, the curves are 
for models with starburst fractions of 100%, 30%, 10%, and 5%. 
The labels along the curves denote the age of the burst in Gyr. 

in the stellar spectra; emission-lines from the interstellar 
medium are not included. While this makes the model 
predictions less meaningful for emission-line galaxies, they 
should be applicable to the E+A galaxies that we are 
primarily interested in. 

Along the track, t he magnitude-weighted stellar age 
ijKauffmann et alJ l2003h from the models agrees qualita- 
tively with the location of the various galaxy subpopula- 
tions. If the E+A population observed near the peak of the 
X-ray emission was created by the equal-mass merger in 
MACSJ0025, these galaxies' approximate age of 0.5-1 Gyr 
(Fig. I17|l would provide a constraint on the time elapsed 
since core passage (Phase C). As interestingly, in the con- 
text of our merger model, is that the regions farther from 
the cluster core along the merger axis are dominated by 
young absorption- line galaxies in which the star formation 
ended about 1 — 2 Gyrs ago. A high fraction of young pas- 
sively evolving galaxies along the merger path is consistent 
with our model in which the transformation of these galaxies 
would have begun in phase A, i.e. at the onset of the cluster 
collision. These time scales are c onsistent with the predic- 
tions from numerical simulations (jBarnesH200lh which show 
that the relative velocity of the two free-falling dark-matter 
halos is significantly reduced by dynamical friction during a 
head-on merger. 

For reference, Fig. [T7] also shows the location of e(a) 
galaxies to demonstrate, qualitatively, that their stellar ages 
derived from this H,5-D n (4000) method are younger than 
those of the E+A galaxies. A more quantitative analysis of 
the e(a) galaxies' locus in this graph would require more 
complex models which are beyond the scope of this paper. 



old galaxy with exponentially decaying SFR (t = IGyr), Solar 
metalicity, and a Salpeter initial mass function. 
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Finally, and not surprisingly, we find the absorption-line 
galaxies in the two sub-cluster core (red symbols in Fig. I17|) 
to be the oldest population, consistent with the expectations 
that the star-formation epoch in elliptical g alaxies at center 
of "p roto" -clusters ends as early as z~ 2(e.g. Blakcslce ct al. 
2006]) . 



6.5 Relation between starburst and cluster 
merger 

Simulations suggest increased sta r-formation activity during 
cluster mergers (e.g. [Bekki 1999). However, most of the ob- 
servational results (e.g. Marcillac et al.l [2007t ) for starburst 
galaxies in clusters are interpreted as being related to infall. 

The spatial distribution of e(a) galaxies in MACSJ0025 
( £15.2.21 and i]5.2.1[) seems to support the infall scenario, in 
that most of the starburst galaxies are located at the cluster 
outskirts. This is consistent with the view that starbursts in 
e(a) galaxies are caused by galaxy-galaxy interactions which 
are highly improbable at the high relative velocities of galax- 
ies closer to the cluster core. 

Nevertheless, evidence of starbursts near the X-ray cen- 
tre is observed in the form of the quenched starbursts in 
E+A galaxies. The physical mechanism that caused these 
bursts is likely to be the same that ultimately quenches 
them, creating the E+A phase of galaxy evolution, namely 
the violent interaction with the dense intracluster medium 
encountered during the merger (b and c, in Fig. ll6[) . The fact 
that no starburst or post-starburst galaxies are found in the 
regions farther from the cluster core, but still well within the 
viral radius, suggests that extreme gas densities are required 
to trigger bursts, which is consistent with the few instances 
in which such dramatic gas-galaxy encounters have actually 
been observed (Cortese et al. 2007; Richard et al. 2010). 
This still leaves us with the question, however, of why no 
active starbursts are presently observed near or within the 
newly formed gaseous core (phase d of Fig. 1 16 p . We argue 
that galaxies from either cluster approaching the core along 
the merger axis, and thus encountering a similarly hostile 
environment as the galaxies at the leading edge of either 
system during first approach, do in fact experience bursts of 
star formation and will, ultimately, also go through an E+A 
phase. Their density, however, and thus the number of such 
ICM-triggered starbursts decreases rapidly as the merger 
proceeds and only galaxies from the increasingly sparsely 
populated outer cluster regions continue to approach. 



7 SUMMARY 

• We present results of a extensive morphological, spec- 
troscopic, and photometric study of the galaxy population of 
MACS J0025.4— 1225, an equal-mass cluster merger, based 
on 212 galaxy spectra obtained with Keck DEIMOS as well 
as colour imaging of the cluster core with HST/ACS and 
wide-field imaging in seven passbands from Subaru Suprime- 
Cam and CFHT/Megacam and CFHT/WIRCAM observa- 
tions. 

• Six E+A galaxies are found in MACSJ0025, all of them 
near the cluster center defined by the X-ray surface bright- 
ness peak, in contrast to other studies which find the E+A 
population to avoid the cluster core. The color of these E+A 



galaxies fall in between those of red-sequence galaxies and 
emission-line galaxies. Compared to E+A galaxies in other 
clusters at intermediate redshift, the luminosity function of 
E+A galaxies in MACSJ0025 appears truncated (no galaxies 
with Mr < —20.7). Their morphology is broadly lenticular, 
i.e., characterized by a dominant bulge and relatively flat 
surface-brightness profiles (Sersic index~ 2). In addition, no 
evidence of tidal features or other asymmetric features is 
found. 

• The fraction of absorption(emission)-line galaxies is a 
smooth increasing (decreasing) function of projected galaxy 
density. A similar correlation is found with distance to the 
cluster center. The fraction of dusty-starburst galaxies (e(a)) 
is dramatically lower inside the virial radius than outside. 

• The global fractions of morphological types (E, S, SO) 
of galaxies in MACSJ0025 are consistent w i th tho se in clus- 
ters at z > 0.5 compiled bv lPoggianti et alj ((2009|), although 
the fraction of SO and S types observed in MACSJ0025 are 
the highest among the sample. While the global SO frac- 
tion is thus not dramatically enhanced, we find the local 
value boosted around the X-ray surface brightness, coincid- 
ing with the excess of E+A galaxies. 

• We propose a simple model in which star-formation ac- 
tivity and morphological transformations are caused by dy- 
namic interactions enhanced and accelerated by the clus- 
ter merger. In this picture, the time since the termination 
of the starburst phase in the observed E+A galaxies natu- 
rally matches the time elapsed since the core passage of the 
merger; about 0.5-1.0 Gyrs ago in the case of MACSJ0025. 
Additional evidence for the causal role of the merger is pro- 
vided by the relatively young age of the absorption-line clus- 
ter members near the leading edges of the two clusters en- 
gaged in the collision. 
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